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1. Introduction 

Our universe appears to be to a high precision fiat and dominated by dark matter (DM) 
and dark energy components of an unknown origin. The case for dark matter is particularly 
strong, with several lines of evidence pointing to its existence from galactic rotation curves 
to CMB and colliding galaxy clusters [1] . The exact value of the present day DM abundance 
to some degree depends on the cosmological model for the accelerated expansion. In 
the current standard model, where the expansion history is attributed to a cosmological 
constant, one finds the best fit values S7 m: o — 0.20 and ~ 0.76 respectively [2]. We shall 
use these values throughout this paper as a reference, although we keep in mind that in 
other models for the acclerated expansion considerably different values of DM density are 
predicted [3]. 

The WIMP (Weakly Interacting Massive Particle) paradigm asserts that DM consists 
of particles which interact very weakly under the Standard Model charges and whose pres- 
ence is felt only through their gravitational interactions. In addition to giving consistent 
explanation to a host of independent cosmological observations, the WIMP paradigm is 
appealing because good WIMP candidates arise naturally in the context of several well 
motivated theories beyond the Standard Model, such as supersymmetry [4] . Another pop- 
ular extension of the Standard Model is Technicolor [5] . In TC the mass patterns of the 
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standard model gauge bosons are explained by new strong dynamics without the need to 
invoke fundamental scalar particles. The masses of elementary fermions are then typi- 
cally explained by the extended technicolor (ETC) interactions [5, 6, 7]. An interesting 
alternative is provided by hybrid models featuring in addition to composite TC scalars 
also fundamental scalar particle which couples to the matter fields via Yukawa interactions 
[8, 9, 10, 11]. 

Simple TC models are severely constrained by the limits on flavor changing neutral 
currents and by the appearance of unwanted additional light pseudo-Goldstone bosons. 
These problems are avoided in walking technicolor theories where the technicolor gauge 
coupling evolves slowly due to a near-conformal behavior. Achieving conformal behav- 
ior using only fundamental representations for technifermions requires a large number of 
techniflavours, and this leads to a large oblique S parameter S ~ 0(1) contrary to the 
observed value S ~ 0. However, considering higher fermion representations increases the 
screening of matter fields and allows approaching conformality with much smaller parti- 
cle content [12, 13]. The phase diagrams for higher representations relevant for walking 
technicolor model building were constructed and candidates for minimal models of walk- 
ing technicolor were proposed in [14, 15]. In particular, it was shown that with just two 
techniflavors in the two-index symmetric (i.e. the adjoint) representation of the SUtc(2) 
gauge group the theory is already close to, or even within, the conformal window; this 
two-color and two-flavor theory is called the Minimal Walking Technicolor (MWT) model. 
In [16, 15, 17] it was shown that MWT is compatible with precision measurements and 
that the MWT composite Higgs boson can be light, with mass of the order of few hun- 
dred GeV [17]. Collider phenomenology of MWT model has been studied in detail e.g. 
in [18, 19, 20, 21, 22], and recently the underlying strong dynamics has been investigated 
on the lattice [23, 24, 25, 26, 27]. 

Despite the simple matter content of MWT model, it possesses rich dynamics and 
implies interesting phenomenological consequences. In particular, the model provides pos- 
sible DM candidates which are the objects of interest for us here. Both technibaryonic [28] 
and leptonic [29] candidates have been investigated earlier. In this paper our goal is to 
update and extend the latter study where the WIMP was identified with a state with the 
quantum numbers of a standard fourth generation neutrino, whose existence is required 
by the internal consistency of the MWT model. However, it was observed already in [30] , 
that such particle interacts too strongly and thereby gives rise to a too small DM-density 
to be consistent with the observations in a universe with an ordinary expansion history. 
In [29] we overcame this problem by proposing that the universe were expanding faster than 
normal during the dark matter freeze-out because of a dynamical dark energy dominance. 
However, the improved constraints from cryogenic dark matter searches have subsequently 
ruled out this scenario [31] . In ref. [29] we only considered pure Majorana and pure Dirac 
neutrinos, while more complicated mass mixing structures are in fact naturally generated 
in the MWT-context [21]. Given the most general mass matrix for the neutrino sector, it is 
easy to arrange that the (lightest) state corresponding to the dark matter has its SM charges 
further suppressed by a small mixing angle 9. This suppression helps in two ways over the 
scanario studied in ref. [29]. First, the reduced annihilation cross section leads to a larger 
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relic density, allowing agreement with the observed DM-density even within the standard 
expansion history of the universe. Second, the reduced couplings lessen the strength of the 
WIMP-nucleon interactions, allowing the model to evade the current cryogenic constraints. 
The stability of our WIMP can be simply imposed by giving a new discrete quantum num- 
ber to all new fields relevant for the DM-sector (or a new ^-symmetry) , in analogy to the 
usual i?-parity in the case of supersymmetric theories. For earlier studies related to very 
heavy neutrinos, see e.g. [32]. 

The idea of a superweakly interacting fourth family neutrino WIMPs was already con- 
sidered in ref. [33] and its relevance for the MWT-scenario was suggested in [29], after 
which the MWT-WIMP scenario including the mixing was considered in ref. [34] . Here we 
improve and extend these analysis in several respects. First, we present a more complete 
TC-implementation of the effective WIMP-scalar interaction sector and couplings. Second, 
we compute the WIMP cross sections much more accurately; where refs. [33, 34] used a 
simple leading term longitudinal WW-approximation for the gauge boson final states, we 
compute the complete annihilation cross sections including all gauge and Higgs boson final 
states. Our final results depend both quantitatively and qualitatively on using the com- 
plete cross sections. Third, we carefully compute constraints on the mass spectrum from 
the latest cryogenic searches and from the oblique corrections implied by the high preci- 
sion electroweak data. Finally we identify an important and to our knowledge previously 
unnoticed effect of relative Majorana phases on the predictions of the model. These phases 
have to be introduced to ensure that the Majorana mass eigenstates are positive definite, 
and to span the entire mass matrix parameter space, one needs to supplement the physical 
masses and mixing angles with the signs of bilinear phase factor products. These products 
enter to WIMP-Higgs and WIMP-Z couplings and then lead to different predictions for the 
relic density and for the oblique parameters as well as to different sensitivity on cryogenic 
searches. In the present case only one nontrivial phase product exists and we give all our 
results for the resulting two independent scenarios. Our final results is that MWT-WIMP 
is a consistent dark matter particle over a range of parameters forming a narrow twisting 
band in the WIMP mass - mixing angle (777.2, 0)-plane, roughly falling within the region 
sin# < 0.4 and 25GeV < m 2 < 45 GeV or m 2 > 80-100 GeV. Most of this parameter 
space will be accessible to the combination of the ongoing or the next generation cryogenic 
experiments and neutrino detectors. In addition, the mass of the heavier neutral state 
mi and the mass of the charged lepton ijie completing the weak lepton doublet, are con- 
strained by the high precision electroweak data such that only narrow finite strips in the 
(777,2, ?7i£;)-plane are allowed for a given cosmologically acceptable (7772, #)-solution. 

The paper is organized as follows: we will introduce the underlying technicolor model 
and discuss the precision electroweak constraints in section 2. In section 3 we describe and 
give our main results of the computation of the relic abundance as a function of the 
mass of the lightest neutrino state, 7772 and the neutrino mixing angle sin 6. Then, in section 
4 we consider the constraints to the model by applying the bounds from the accelerators, 
cryogenic searches and from neutrino detectors such as Super-Kamiokande and IceCube, 
including estimates for the future reach of these experiments. Our conclusions and outlook 
are presented in section 5. Finally, in the appendix we present the relevant details of the 
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cross section calculations. 



2. Minimal Walking Technicolor and the Fourth Generation of Leptons 
2.1 Model Lagrangian and mass terms 

In the MWT model to be studied here, the electroweak symmetry breaking is driven by 
the gauge dynamics of two Dirac fermions in the adjoint representation of SUtc(2) gauge 
theory. The key feature of this model is that it is (quasi) conformal with just one doublet 
of technifermions [14]. However, since the technicolor representation is three dimensional, 
the number of weak doublets is odd and hence anomalous [35] . A simple way to cure this 
anomaly is to introduce one new weak doublet, singlet under technicolor and QCD color 
[14, 15] in order not to spoil the walking behavior and to keep the contributions to the 
oblique corrections as small as possible. Hence, the model requires the existence of a fourth 
generation of leptons. The anomaly free hypercharge assignments for the new degrees of 
freedom have been presented in detail in [15]. Here we simply note that there exists an 
assignment which makes the techniquarks and the new lepton doublet appear exactly as 
a regular standard model family from the weak interactions point of view and it is this 
assignment that we will consider throughout this work. 

The technicolor sector of the theory confines at electroweak scale and is described 
better in terms of a chiral effective theory than using the fundamental techniquark and 
-gluon degrees of freedom. The global symmetry breaking pattern is SU(4) — > SO(4), with 
nine goldstone bosons. Three of these are absorbed into the longitudinal degrees of freedom 
of the weak gauge bosons, and the low energy spectrum is expected to contain six quasi 
Goldstone bosons which receive mass through extended technicolor interactions [5, 36, 37]. 
Their phenomenology has been investigated elsewhere [18, 19]. Here we will use the result 
obtained in [18] that the conservation of the hypercharge allows only the effective SM- 
like scalar Higgs to couple to the fermions. Hence, we re-introduce a (composite) Higgs 
doublet into the theory. We denote the left-handed fourth generation lepton doublet by 
L = (Nl,El) and the right-handed SUx,(2) singlets as Er and Nr. To take into account 
the effects of the scalar sector on these leptons up to and including dimension five operators, 
we consider following effective interactions 



where H = ir 2 H* and A is a suppression factor related to the more complete (ETC) 
ultraviolet theory. The first terms in Eq. (2.1) lead to the usual (Dirac) mass for the 
charged fourth generation lepton, and the remaining terms allow for more general mass 
structure of the fourth neutrino. After symmetry breaking the effective Lagrangian (2.1) 
gives rise to a neutrino mass term: 



4iass = (yL L HE R + h.c.) + C D L L HN R 

+ ^(L c H)(H T L) + ^(H^H)N R N R + h.c. 



(2.1) 




(2.2) 
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where ul = (Nl, N r c ) t , txid = Cl>v/V2 and Ml,r = Cl,rv 2 /2A, where tj is the vacuum 
expectation value of the effective Higgs field. The special cases are a pure Dirac and a pure 
Majorana neutrino which are obtained, respectively, by discarding dimension five operators 
and by removing the right handed field Nr. The most general mass matrix contains, even 
after the field redefinitions, one complex phase. However, in this paper we shall restrict 
ourselves to the case of real mass matrix. Mass eigenstates are two Majorana neutrinos 
which are related to the gauge eigenstates by a transformation 

N = On L + pO T n c L , (2.3) 

where N = (N u N 2 ) T and O is an orthogonal 2x2 rotation matrix, where the associated 
mixing angle is 

tan2fl= J m \ f . (2.4) 
Mr- M L v ' 

The phase-rotation matrix p = diag(pi,p2) is included above to ensure that the physical 
masses mi 2 are positive definite. Indeed, the eigenvalues of the mass matrix in (2.2) are 

A± = 1 -(m l + Mr ±^{M l - MrY+AtuI ) . (2.5) 

Because the signs and relative magnitudes of Ml,r and mu are arbitrary, the eigenvalues A± 
can be either positive or negative. However, choosing independent phases as p± = sgn(A±) 
we get positive m± = \X±\ as required. For our purposes it will be convenient to express 
everything in terms of the physical mass eigenvalues 777 1 > 777,2 and the mixing angle sin 9 
instead of the Lagrangian parameters Ml, Mr and mo- While working with physical 
parameters has obvious advantages but the downside is that the connection between the 
physical and the Lagrangian parameters is not always straightforward. 

2.2 Accounting for the mixing phases 

Let us now discuss in detail the role of the phase-rotation matrix p introduced in the 
previous subsection. There are two distinct but equivalent ways to treat these phases. 
First, one may embed them into the definition of the mass eigenstates which leads to 
unitary rotation matrix in Eq. (2.3). Then phases appear explicitly in the interaction 
vertices and related Feynman rules while the Majorana eigenstates satisfy n c = n, making 
the field operator /9-independent. This formulation was considered for this model in [21]. 
The second posibility, which we shall adopt here, is to retain the definition (2.3) with 
an orthogonal rotation. This has the consequence that the Majorana eigenstates satisfy 
7i c = pn i.e. the phases now appear explicitly in the field operators corresponding to 
the Majorana eigenstates. Then the phases have to be properly taken into account in 
the contractions and propagators when evaluating various processes. This issue will be 
considered to some detail in the appendix. 

Let us now describe the structure of the parameter space and the role of the phases. 
In Fig. 1 we show the parameter space along a plane of some nonzero value of mo. Let us 
first suppose that Mr + Ml > 0. Then the eigenvalue A+ is (2.5) is always positive and 
larger than |A_|. We choose to denote by N± the heavier state, so that p\ = p + = 1 and 
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Figure 1: Shown is a constant-™/) slice of the mass parameter space, which is divided into 
physically disjoint regions by the sign of the effective phase p+p-. The physical parameters are 
invariant under the mirroring symmetry with respect to the line Mr = —Ml- 



our WIMP is always associated with N2 with the mass m<i = m_. The sign of P2 = p~ is 
defined by the ratio of m 2 D and MlMr: p2 = — 1 if m 2 D > MlMr, and otherwise P2 = +1. 
Individual phases are not observable, but relative phase corresponding to the product of 
the two is: 

P1P2 = P+P- = P12 • (2.6) 

In particular this phase will show up explicitly in the various couplings of the mass eigen- 
states. The situation with respect to the product phase is again illustrated in Fig. 1: p+p- 
is positive in the upper right corner separated by the black solid line and negative else- 
where in the region Mr + Ml > 0. From Eq. (2.4) we see that tan 29 becomes infinite 
along the the line Ml = Mr, corresponding to a maximal mixing \9\ = tt/4. We define 
the mixing angle to be zero in the limit Mr — Ml >> |m£>|, whereby < \9\ < tt/4 below 
and 7r/4 < \9\ < tt/2 above the line of maximal mixing. The sign of the mixing angle is 
determined by the sign of mo- Physically, in the region below the line of maximal mixing, 
our WIMP is predominantly a sterile right chiral state while in the region above the line it 
is predominantly a left chiral state with ordinary weak interaction strength. It is thus clear 
that the area of most interest for us is the rightmost quadrant bounded by the diagonal 
lines in the phase space in Fig. 1. 

Now consider the case Mr + Ml < 0. In this area we have that A_ is always negative 
and I A_ I > |A+|. Thus, in this region we have to associate our heavier state (which we 
always choose to label with N\) with A^_, while the WIMP is always the now lighter 
A^2 = N + . Here we have always p\ = p_ = — 1 while the sign p\ = p + again depends on the 
relative magnitude of MlMr and m 2 D ] the resulting division to distinct areas according to 
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the sign of pu = p-p+ is again shown in Fig. 1. There is an obvious symmetry in the 
phase space about the reflection along the line Ml + Mr = 0. Indeed, all that happens in 
this reflection is that the eigenvalues A± and their associated eigenstates exchange roles, 
but all physical parameters 9, m±, m 2 and p\ 2 remain invariant. That is, the reflection 
corresponds to a mere relabeling 1 <-> 2 as the two regions can be mapped to each others 
by a redefinition of the phases of the states. We shall thus concentrate on the region 
Mr + Ml > from now on. Nevertheless, for each triplet of physical mass and mixing 
parameters mi, m 2 and 6 the original parameter space contains two physically distinct 
solutions labelled by the relative phase pn. In what follows, we shall always give the 
results for both possibilities. 

We conclude this subsection by noting that all the typical special cases are contained 
within the pn = — 1 portion of the parameter space: the pure left- or right-handed Ma- 
jorana states correspond to Mr = and Ml = axes in the plane mo = 0, respectively, 
while the Ml = and Mr = planes for nonzero mo correspond to usual seesaw sce- 
narios. Finally the Dirac limit corresponding to the axis Ml = Mr = is also contained 
only in the pu = — 1 domain. It is perhaps due to this reason that the other domain with 
pn = +1 has so far gone unnoticed in the literature. 

2.3 Couplings 

For the analysis of the relic density we need the couplings of the neutrino mass eigenstates 
to the weak gauge bosons and to the Higgs boson. These are easily found out by applying 
the appropriate phase- and rotation transformations defined in the previous section. We 
shall write down only the terms relevant for our calculations. For the Z and W ± bosons 
we find that 

W+N L -fE L = sin^ W+Nil^El + ■■■ 
Z^Nl^Nl = sin 2 9 Z^Nil^Nil 

+ \ sin 29 Z^ (N 1l ^N 2 l + N 2L ^N 1L ) + • • • , (2.7) 

where the omitted terms contain interactions of the heavy N\ field only. These couplings 
are diagonal in the mixing and therefore do not involve the phase factor pu. However, 
neutral current involves mixing and these couplings do depend on pu. One finds: 

Ni-fZpPLN! + N^Z^PlNi = Ni((3 + cr^-fZ^, (2.8) 

where 

a=|(l + Pi 2 ) and (5 = \{l - p 12 ) . (2.9) 

Thus, for pu = — 1 the neutral current interaction of our WIMP is purely axial vector 
and for pu = +1 purely vector. Usually in the literature dealing with the interactions of 
Majorana neutrinos, only the first possibility is mentioned. 

The effective interaction terms involving the Higgs and the lighter neutrino eigenstate 

are 

Cnh = ( C^hNiNi + C^hN^a - ^)N 2 
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+ C £h 2 N 2 N 2 ) + n ^-h 3 + -- - , (2.10) 

where we have again omitted the interaction terms which do not contain N 2 and hence are 
not needed in our analysis. The factors a and (3 are defined in Eq. (2.9) and the factors 
C 22 , C21 and C22 are defined in the left panel of Table 1. 





Scenario I 


Scenario II 




1 - \ sin 2 29 R- 


sin 2 # 


U 21 


-\p 12 sm46 R- 


\p\ 2 sm29 R+ 


°22 


\-\ sin 2 29 R_ 


\ sin 2 9(1 - cos 2 OR-) 



So far our model building has 
rested on the assumption that the 
(composite) Higgs is the sole source 
of elementary fermion masses. How- 
ever, this is not necessarily the case, 
and there may be other (composite 



Table 1: Coefficients of the Lagrangian (2.10) for two the or even fundamental) scalars whose 
distinct mass generating scenarios described by Eqs. (2.1) condensation i eads to mass terms 
and (2.11). We have defined R ± ee 1 ± p 12 ^. ^ ^ ^ To mustrate 

such possibilities, we consider as an 
alternative to the model Lagrangian (2.1), the case where the right-handed neutrino mass 
originates from a Standard Model singlet scalar field S. 

£lLs = {vLlHEr + h.c.) + C D L L HN R 

+ ^(Z C #)(# T L) + CrSNrNr + h.c. (2.11) 

This model is similar to the usual see-saw neutrino mass generation mechanism, although 
here the singlet S does not need to be a fundamental scalar. To specify the model com- 
pletely one should give a potential for S. However, none of the parameters of this potential 
are needed in our analysis; we may assume that the vacuum expectation value for S is 
generated through interactions with the Higgs, i.e. we do not need additional sources of 
spontaneous symmetry breaking for the dynamics of the 5-field. In what follows, we will 
refer to the scenario with just only the doublet Higgs field as Scenario I and to the case 
with Higgs and a singlet scalar as Scenario II. The interactions between the Higgs and the 
neutrino can be generically described by the Lagrangian (2.10) for both scenarios. The 
precise form of the coefficients is given in Table 1 . 

Let us conclude this section by noting that to keep the WIMP stable, we have to 
exclude all renormalizable operators of the form 

W+N tLl n R , Z^l^Nr, L l H£ r , ... (2.12) 

which would mix the new WIMP-sector to the other standard model fields L^, £r, or v^r 
where t = e, p, r. The simplest way to achieve this is to assume that all new fields relevant 
for the DM-sector, Ll, Er and Nr share a new discrete conserved quantum number. 
Then all allowed Lagrangian terms must obey a new ^-symmetry under the exchange 
of these new fields. This construction is analogous to the i?-parity imposed in Minimal 
Supersymmetric Standard Model to guarantee stability of protons (and consequently of 
the LSP). 
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2.4 Oblique constraints 

The fourth generation of leptons is constrained by current accelerator data. From LEP 
we know that the charged lepton E has to be more massive than the Z boson and if 
the fourth generation neutrino has standard model interaction strength, it needs to be 
heavier than Mz/2 in order to evade the constraint from Z-pole observables. In the case 
of neutrino mixing considered in this work, the lighter state can have a substantial right- 
handed component and hence interact only very weakly. This, as we shall see, can allow 
this state to escape the LEP bounds even when its mass is less than Mz/2. In addition 
to these direct bounds, the parameters of the fourth generation leptons are constrained 
by oblique corrections, i.e. due to their contribution to the vacuum polarizations of the 
electroweak gauge bosons. These contributions are conveniently represented by the S and 
T parameters [38] . 

The oblique corrections in MWT model with the general mass and mixing patterns 
considered here have been studied in detail in [21]. We also note that there exists two 
extensive fits performed by the LEP Electroweak Working Group (LEPEWWG) [39] and 
independently by the PDG [40]. Both fits find that the SM, defined to lie at (S, T) = (0, 0) 
with m t = 170.9 GeV and ran = 117 GeV, is within la of the central value of the fit. The 
two fits disagree on the central best-fit value: LEPEWWG finds a central value (S, T) = 
(0.04, 0.08) while including the low energy data the PDG finds (S, T) = (-0.04, 0.02). Since 
the actual level of coincidence inferred from these fits depends on the precise nature of the 
fit, we allow a broader range of S and T values, roughly corresponding to the 3a contour. 
Concretely we require < S < 0.2 and < T < 0.5. From the results of [21] it can be 
inferred that these values can be accommodated easily within the parameter space of the 
leptonic sector. 

In this work we will supplement the constraints on the parameter space by requiring 
in addition to the saturation of the precision constraints that the lightest neutrino mass 
eigenstate provides the correct relic density to match the observed DM abundance. Hence, 
we adopt the following strategy: We will first determine the mass of the lightest eigenstate, 
777,2 and the corresponding mixing angle sin 9 such that the relic abundance f2jv 2 (777,2, sin 9) ~ 
0.2. This constraint between 7712 and sin# will be practically independent of the values of 
the remaining mass parameters m\ and txie in the leptonic sector provided 7772 <C mi,mE- 
Then, for each 7772 and sin 9 constrained through £In 2 we will determine the allowed values 
of mi and 777£ by requiring that resulting contribution to S and T is within the bounds 
quoted above. We will therefore now move to describe the evaluation of the relic density 
and the constraints from earth-based direct dark matter searches and return to the oblique 
corrections and constraints in Sec. 4.3. 



3. Relic Density 

The relic abundance S7tv 2 is computed in the standard way. We start from the Lee- Weinberg 
equation for the scaled WIMP number density [41]: 

df(x) _ (va)mlx 2 2 2 



dx H 



(f(x)-fL(x)), (3-1) 
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where m 2 is the WIMP mass and we have introduced the variables 

( \ 1//3 
/We^, and x = — ^ — , (3.2) 
•SB m 2 

where se(T) is the thermal entropy density at the temperature T. Given se, the Hub- 
ble parameter H(T) = (87rp(T)/3M| 1 ) 1 / 2 and the average WIMP annihilation rate (va), 
Eq. (3.1) is easily solved numerically. Here we assume that the expansion of the universe 
follows the standard adiabatic expansion law, so that H and se can be computed from 
their standard thermal integral expressions. Typically our WIMPs are freezing out at 
T ~ 0(1 — 10) GeV, so that the uncertainties in se arising from the QCD phase transition 
play no role for us. After the present ratio of ^-number-density to the entropy density /(0) 
is found from Eqn. (3.1), the fractional density parameter Qn 2 of the Majorana WIMPs 
becomes 

n N2 ~ 5.04 x 10 5 m 2 /(0) . (3.3) 

From Eq. (3.1) one sees that the relic density /(0) depends essentially on the ratio (va)/H. 
The smaller this quantity is, the less time the WIMPs can remain in thermal equilibrium 
and thus the larger is their relic abundance. One can show (see e.g. [42]) that the depen- 
dence is in fact almost linear: Vt^ 2 ~ H/ (va). With the standard expansion history of the 
universe H is fixed, and so the characteristics of the solution are entirely dictated by the 
average cross section (va). 

We compute (va) in the Maxwell-Boltzmann approximation [43]: 

1 f°° - fs 

where Ki(y)s are modified Bessel functions of the second kind and s is the usual Mandelstam 
invariant. This approximation for the collision integral is accurate to within a few per 
cents for massive neutrinos [44]. For the total cross section <r tot we considered the A^A^ 
annihilation to the final states including all open fermion, gauge boson and scalar channels 

N 2 N 2 -► //, W + W~, ZZ, ZH° and H°H° . (3.5) 

We omitted annihilations to technifermions because these rates would be just a small 
correction to already subleading fermionic channel. Above, H° is the effective, light "SM- 
like" Higgs state appearing in the mass operators (2.1) and (2.11). We did not include the 
SS final states in the scenario II, assuming that the new scalar S is heavy. We computed the 
complete cross sections for each channel shown in (3.5) without further approximations and 
performed all s-integrals numerically. The iV2-gauge boson and A^-Higgs couplings needed 
in these computations were given in section 2.3. In this work we make the assumption 
that the heavier, unstable neutrino N± has already decayed and is no longer present during 
the N 2 freeze-out. That is, we assume that the particle spectrum during the freeze-out is 
just the usual Standard Model particle spectrum and the annihilating WIMP. More details 
about the computation of the cross sections are given in the appendix. 
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It is clear that the most important parameters setting the scale of {va), and hence 
that of ^tv 2 , are the WIMP mass mi and the mixing angle 9. We have therefore displayed 
our main results as £In 2 contours in the (777.2, sin 0)-plane. The results are only very weakly 
dependent on the sign of the mixing angle however, and we will always fix sin 9 positive in 
what follows. The charged lepton E appears only as a virtual state in the i-channel W + W~ 
process and so our results depend only very weakly onm^ as well. The dependence on m\ 
can in principle be strong 1 , but it turns out that in the region of interest, where the relic 
density can be large enough, and the observational constraints are satisfied (at relatively 
small mixing angles), also the dependence on m\ is weak. For definiteness we have set 
the = 777i = 27772- On the other hand, our results are very sensitive on the mass of the light 
composite Higgs particle 777/7 as we shall see below. In addition to these mass parameters 
and mixing our results depend on the relative phase factor pw This dependence can be 
seen explicitly in the couplings derived in section 2.3 and in the matrix elements given in 
the appendix. Finally, we have chosen to consider two different mass generation schemes 
in this paper; the first one (scenario I) using only the light composite H° and the other 
(scenario II), where the Majorana mass of the right-chiral state is generated by a light 
singlet S. The AVHiggs couplings in particular depend very sensitively on the choice of 
the scenario. 

To summarize: Our model predictions for the relic density Oat 2 are most dependent 
on parameters 7772 and sin 9. They are essentially sensitive to ran, pyi and on the mass 
generation scenario, while they have only a small subleading dependence on the masses m\ 
and 777£. 

We show our results 2 in Fig. 2 where we used mu = 200 GeV and in Fig. 3 where 
we took run = 500 GeV. In both figures the left panel corresponds to the relative phase 
P12 = +1 and the right pane to p\2 = — 1. In all figures we show the contours fijv 2 = 0-19 
and £In 2 = 0.23 corresponding to the boundaries of the region where Q,n 2 1S consistent 
the WMAP [2]. The thick dashed lines mark the contours for the scenario I and thick 
solid lines for scenario II. In addition we have shown by thin dashed lines the contours 
corresponding to Hjv 2 = 0.1 and fijv 2 = 0.5 in the scenario II. We have also shown the 
regions of the parameters that are already excluded by the direct LEP-constraints (light 
yellow area) and by the XENON10 cryogenic dark matter search (dark red area), which 
are currently the strongest constraints on the model. We have also shown (dash-dotted 
curves) the sensitivity of the future XENON100 experiments on the model parameters. The 
observational constraints will be discussed in more detail in the next section. We finish 
this section by building an intuitive understanding of the relic density contours shown in 
Figs. 2-3. 

Let us first concentrate on region 7772 < Myy , where the dominant annihilation channel 
is the one to standard model fermions. This cross section is simple enough to be given here 



x Our generic system includes also the pure Dirac limit, where the predictions are of course entirely 
different from the ones shown here. 

2 These figures are analogous to Figs. 1-2 of ref. [33], with the identification g' /g = sin#. Note the typos 
in the powers of g'/g in Eqs. (5-7) of [33]: they should be ±4 instead of ±2. 
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Figure 2: Shown are constant Ojy 2 contours as a function of mass and mixing angle. The area 
between the contour lines marked by 0.19 and 0.24 is consistent with the WMAP results for the 
dark matter density parameter. In both panels m# = 200 GeV, thick solid lines correspond to 
the scenario II and thick dashed lines to mass scenario I. In left panel we took p\2 = +1 and in 
right panel p±2 = — 1. Thin dotted contours show additional contours for the scenario II. Yellow 
(light shaded) area is excluded by the LEP limits and the red (dark shaded) area is excluded by 
the XENON10 limits. The red dash-dotted lines show the predicted sensitivity of the XENON100 
(upper line) and XENON100 update (lower line) as given in ref. [45]. 



explicitly: 

G 2 F m^ P f f f sin 



'M^ffi') = -teT-pfb (^H^l 9f(s,m 2 ,rn f ) 

+M(C% 2 ) 2 \D H \ 2 s 2 f3 2 p 2 f m 2 m 2 } . (3.6) 

Here 0w is the Weinberg angle, s = s/rriyy, m 2 . * = m| ^/m 2 v , fix = (1 — Am 2 x /s) l l 2 and 
\t) x \ 2 = m^\D x \ 2 for X = Z,H, with 

Dx = 2 ~^ , (3.7) 

s — m x + ii xf^x 

The factor C 22 accounts for the different coupling strengths of the lightest neutral particle 
to the SM-like Higgs as indicated in Table. 1. Finally N^, = 1 for leptons and N^, = 3 for 
quarks and 

gf(s, mi, rhf) = (§s 2 — |sm| + 2m 2 m 2 f S J (v 2 + a 2 ) + 4mj(s — Qm 2 )(v 2 — a 2 ) , (3.8) 

where vj = — 2Qfs'm 2 9w and aj = T^t where T3 is the isospin and Q is the charge of 
the fermion. Because £In 2 is inversely proportional to the cross section, the predicted relic 
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Figure 3: Same as in figure 3, but now for mjj = 500 GeV. 

density is smallest near the Z-pole for a constant mixing angle. Conversely, to keep Qjv 2 
constant, we have to compensate for the Z-pole by decreasing sin#. This tradeoff results 
in the characteristic dip at m 2 ~ Mz/2 apparent in Figs. 2-3. The Higgs contribution, 
proportional to |-D//| 2 in Eq. (3.6) is always subdominant to the Z-exchange, because of 
the small Yukawa couplings of the light SM-fermions. This explains why our results for 
different scenarios and for different choices of phase factor p\ 2 are almost identical for 
m 2 <, M w . 

For M\y < m 2 ^ Mfj the total cross section <7 t ot is dominated by the annihilation 
into gauge boson final states. All channels are important quantitatively, but the most 
important qualitative features of the results are set by the Z- and ff-mediated s-channel 
annihilations to longitudinal charged gauge bosons. The cross section for the former of 
these sub-processes near the mass threshold is 

a(N 2 N 2 ^Z*^ WfWT) « ^Isin^fl - ^) . (3.9) 

6-7T V S J 

If this was the only contribution, we would get (see e.g. [33].) (va) ~ G^m^Tsin 4 9/2tt, 
and eventually 

TeV 2 

« 1.2 x 10~ 3 2 . (3.10) 
sm 

This predicts a slope sin# ~ 0.3(TeV/m) _1 / 2 for a constant J7jv 2 ~ 0.2; the actual slope 
has a smaller coefficient because of the contribution from ZZ and ZH final states. This 
scaling was observed both in [33] and in [34] and combined with the XENON10 bound it 
would exclude masses m 2 5^ 300 GeV. However, this simple scaling is completely changed 
when one includes interactions with a relatively light Higgs field. The most striking effect 
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occurs with the scenario I where no solutions are found when mi > Myy. The is due to the 
Higgs mediated s-channel annihilation to the W /+ VF~-final state. This rate is proportional 
to the Higgs- WIMP coupling which in scenario I is not suppressed by the mixing angle. 
The complicated behaviour {e.g. with the curves bending back towards smaller mass with 
decreasing angle) in Figs. 2-3 are caused by the complex phase- angular dependence of C 22 
and by the m^-suppression in the Higgs propagator. In scenario II all Higgs couplings 
are proportional to the mixing and so the Higgs effects are more moderate. Nevertheless, 
for a relatively light Higgs mass, the above mentioned s-channel interaction creates the 
characteristic H-pe&k at m 2 ~ mff/2 similar to the Z-peak seen at m 2 ~ mz/2. Beyond 
the Higgs peak the s-channel exchange becomes less important and one returns to the 
s'mO ~ m _1//2 -slope. For m 2 ^ ran the behaviour changes again (only the scenario II is 
interesting here), because of the opening of the HH- final states. In particular the four- 
point contact term has a large coupling, partly because of combinatoric factors and partly 
because of the mass dependence C 22 ~ mi/m2 at small mixing. 

To summarize, we find that the superweakly interacting MWT-WIMP is a good dark 
matter candidate for a range of mass and mixing parameters. The preferred scenario is 
the one where the light neutral mass comes from a new singlet field. The WIMP of the 
scenario I, which uses only the field H for the mass generation, is restricted to have the 
mass 30GeV <> m 2 <, 50 GeV or m 2 ~ M w + 0.15(m H - 200) GeV, and sin6» £ 0.4. 
In the scenario II employing a new singlet S, the high mass solution is opened up to a 
region m 2 <: 100 GeV and sin <, 0.4 — 0.5 with exact bounds depending on the mass 
of the H-field. Let us note that the coincidence regions where the WIMP mass is close 
to some resonance, m 2 ~ Mj/2, are generic sweet spots for the superweakly interacting 
dark matter scenarios; in these cases DM particles can have a sufficiently large annihilation 
rate despite a very small coupling, which makes them difficult to detect in direct searches. 
Nevertheless, the near future XENON100 updates will be able to either detect, or all but 
rule out the WIMP proposed here, as shown by the predicted sensitivity of these searches 
in Figs. 2-3. 

Let conclude this section on a note on the MWT-WIMP scenario without mixing 
considered in [29] . There a sufficiently small relic density was not arranged by a superweak 
cross section like here, but by a large expansion rate of the universe in the context of 
quintessence-like dynamical dark energy model. As we have mentioned, this scenario is 
essentially ruled out by the XENON10 data. This is in fact evident from our Figs. 2-3, 
where the line sin 9 = 1 would correspond to any desired O^-contour in the pure Majorana 
case of ref. [29]. The Dirac case is even more constrained. If we changed the expansion 
history of the universe in the present mixing scenario, the constant f^-curves would tend 
to move towards a larger mixing, and hence closer to or within the excluded region. In 
other words, changing the expansion history would change the predictions in a way that 
would further restrict the allowed parameter space. 

4. Constraints 

In previous section we presented our main results, including the presently strongest ob- 
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servational constraints. Here we will analyze in more detail how these constraints on the 
(mjv, ^-parameter space arise. Moreover, we will consider the present sensitivity and the 
future reach of the other dark matter searches, such as CDMS [46], Super-Kamiokande [47] 
and IceCube [48] experiments on the parameters of our model. Finally, we will also com- 
pute the constraints from the precision electroweak data. Here the strategy is somewhat 
different, because the oblique corrections are sensitive mostly to the more strongly inter- 
acting (heavier) neutrino and to the new charged lepton. Fortunately the relic density 
analysis depends on these parameters only very weakly, as we discussed in the previous 
section. Oblique corrections will therefore be computed in the (mi, m,£)-plane for a number 
of cosmologically acceptable (mjv, #)-pairs. 

4.1 LEP-limits 

Let us start by reviewing the constraint coming from the LEP-measurement of the Z-boson 
decay width to an invisible sector [49] . Written in terms of the number of massless light 
neutrino species the allowed width to the invisible sector [50]: 

N v = r ( Z -" mv ") = 3.00 ± 0.08. (4.1) 

1 theoryl^ ~ ► VV) 

The known three light neutrino families almost saturate this bound, leaving only the de- 
viation of the experimental value N„ from three for the new invisible particles. At the 
Z-peak, the relation between the Z-boson decay to a massless SM-neutrino and to our 
heavy neutrino is T(Z — > N2N2) = sin 4 ^/?! x T(Z — > vi>), so that from Eqn. (4.1) we get 
the constraint: 

0.08 > N v - 3 = sin 4 9 x /3g , (4.2) 

where is the neutrino velocity factor with s = m 2 z . The exclusion contours corresponding 
to this limit were shown by the yellow shaded area in the Figs. 2 and 3. 

4.2 Cryogenic and other indirect detection limits 

Adapting the constraints from a given dark matter search to a particular WIMP model is 
not always straightforward. The problem is that the WIMP-nucleus interactions depend 
on the precise form of the WIMP-nucleus and of the WIMP-nucleon couplings as well as on 
the proton and neutron structure functions within the target nucleus. In an ideal situation 
one would compute precise experiment-specific expected count rate for the WIMP model 
one is interested in: 

proc. isot. 

N = m d rJ2^Z^a (4.3) 

a i 

proc. isot. y fE max pv csc 

= — E E ^ i *> *»« i UQ) *-AW^(9Wfl) • 

This is a fairly complicated expression with many different parameters: m^ is the mass 
of the detector, r the detector exposure time, R4 is the fractional count rate, Yi is the 
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Factor 


Spin-dependent (axial vector) 
coupling 


Spin-independent (scalar) 
coupling 


a 




4 2 


F?{Q) 


S(q) 
5(0) 


[^f)]%xp[-(„) 2 ] 



Table 2: Shown are the spin-dependent and spin-independent neutrino- nuclide coupling factors. In 
spin-dependent case the factors Cj and Ff(Q) are from [31] the equations (2) and (3) respectively. 
In spin-independent case Ff(Q) is from [52]. (q = ^/2Qmjv i -)) 

fractional abundance and rn/V; is mass of the target nuclide of a specific isotope I and 
P2 = 0.3GeV/cm 2 and mi are the local WIMP density and mass. Q is the energy deposited 
on the nuclide in the elastic collision, Et is the detector threshold and £ max is the maximum 
energy which nuclide can get from the impact. Furthermore, v is the WIMP velocity, 
fi(v) is the WIMP velocity distribution in the earth frame [51, 29], v^ n {Q) = 
corresponds to the minimum velocity of the WIMP, v csc is its escape velocity from the 
galaxy gravitational potential and /Xj = mivn^Jimi + mjyj is the reduced mass of the 
WIMP and the target nuclide. Finally, the cross section Oi a is given by 

Via = Ciatfaoa (4.4) 

where Cj a is an enhancement factor shown in table 2, and ao a is a dark matter model 
dependent constant which is independent of the isotope. The formula (4.3) involves also 
the nuclear isotope, interaction process and nuclear model specific functions Ff a {Q) and 
the process and experimental setup dependent function Xa{Q) that characterizes for ex- 
ample the detection efficiencies and cuts used to reduce the background. Unfortunately 
these functions are often insufficiently detailed by the experiments, making it difficult or 
impossible to use Eq. (4.3) directly. Conversely, it is not feasible for experiments to put 
forth specific constraints on all different DM-models and so, in a compromise between 
accuracy and generality, the observational limits are typically expressed in terms of pure 
WIMP- proton and WIMP-neutron cross sections. 

The most important quantity affecting the constraints comes from the nuclear spin 
dependence of the WIMP-nucleon interaction. Interactions which require a nucleon spin- 
flip to proceed, are called spin-dependent interactions. Our WIMPs are Majorana particles 
with an axial vector coupling to Z and a scalar coupling to H. Thus, the Z-mediated 
interaction process for our model is spin-dependent, while the //-mediated process is a 
spin-independent one, which does not require a nucleon spin-flip to proceed. All WIMP- 
nucleon interactions are very soft and so the spin-dependent interactions are only seen 
with nuclei with (an) unpaired valence nucleon(s). Because WIMPs couple differently to 
protons and neutrons, the identity of these valence nucleons (whether mostly a proton or 
a neutron) gives rise to a further sensitivity on nuclear physics details. At any rate, to get 
the most accurate bound on a given model one should sum the spin-independent and spin 
dependent processes incoherently (because these processes have different final states) in the 
predicted count rate as indicated in Eq. (4.3). However, in practice one typically assumes 
that WIMP has either purely spin-dependent or purely spin-independent interactions. 
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For our model, the best current limit comes from the cryogenic dark matter search 
XENON10 experiment [31]. Fortunately, this experiment has given their (spin-dependent) 
constraints not only for the standard proton and neutron interactions, but also for a stan- 
dard model 4 th family Majorana neutrino. Better yet, reference [31] explicitly plots the 
expected count rate in their detector for this case. Now, the spin dependent count rate N 
predicted for our WIMP through Eq. (4.3), differs from the standard model case only by 
a simple scaling of the cross section factor ao in Eq. (4.4): <5"o,SM — ► <5"o,Mix = sin 4 ^cj^sm- 
Using this information we can convert the XENON10 results for a 4 th family SM-neutrino 
to an upper limit on the mixing angle as a function of mass: 



The function iV~sM(77i2) was read from the left panel in Fig. 2 of ref. [31] and the function 
Mimit( m 2) was approximated by a linear interpolation between the values of NsM.('m2) 
at the high and low mass ends of the SM-exclusion region in the same figure. The red 
areas in the Figs. 2 and 3 show the excluded region corresponding to the upper limit (4.5). 
Finally, the future reach of the XENON100 updates can be simply estimated by scaling 
the above limiting angle (4.5) by the fourth root of the mass-day exposure ratios of the 
present and future searches: sin#i im — ► sin OumiExeio/Eupg) 1 ^. We use Exeio = 136 kg- 
days for XenonlO, Exeioo = 6000 kg-days for XenonlOO and ExeWOu = 60000 kg-days for 
the XenonlOO upgrade [45]. The results are displayed by the red dash-dotted curves in 
Figs. 2 and 3. 

As mentioned above, our WIMP has also spin-independent interactions mediated by 
the Higgs field, which were not accounted for in the treatment leading to the constraint 
(4.5). Let us check that the limit coming from this process is really subdominant. This 
is not immediately clear, because the spin- independent channel is coherent, and it is thus 
enhanced with respect to a spin-dependent one by a factor ~ A 2 F 2 (Q) where A is the mass 
number of the target nucleus and F 2 (Q) is a shadowing correction shown in table 2; for 
the XENON10 experiment the enhancement is roughly ~ 10 4 . The spin-independent limit 
from XENON10 is given in [53] in terms of the WIMP-nucleon cross section. In our model 
this cross section (in the zero momentum transfer limit 3 ) is given by 



where n refers to a nucleon and fi n is the WIMP-nucleon reduced mass. Note that without 
the C^-factor the cross section (4.6) would coincide with the SM Higgs-neutrino scattering 
cross section. Our results are shown in Fig. 4. The red dashed line corresponds to the lower 
limit for the WIMP-nucleon cross section from XENON10, and the blue dashed line shows 
the same constraint from the CDMS experiment [46]. The solid line shows WIMP-nucleon 
cross-section of Eq. (4.6) corresponding to the parameters mi and sin #(7712), constrained 

3 The zero momentum transfer approximation could be improved by the methods of ref. [55]. These 
corrections are quite small however, and they would not affect the final conclusion. 




(4.5) 




G 2 F nl m\m\ 



(4.6) 
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Figure 4: Shown are spin-independent WIMP-nucleon cross sections. Dash dotted blue line shows 
the latest CDMS constraint [46] and the red dashed line corresponds to XENON10 limit from [53]. 
The solid black lines show the predicted cross section yielding f2jv 2 = 0.214 in scenario II with 
ma = 200 GeV and p\2 = +1 (left panel) and with mjj = 500 GeV and pi2 = — 1 (right panel). 
The dotted gray line shows the standard, mixing angle independent S1jv 2 = 0.214-cross section 
for comparison. Finally, the light dashed (green) lines show the sensitivity of the XENON100 
experiments and the light dotted (green) line that of the XENON1T experiment [45]. For the 
production of the experimental curves here and in fig. 5 we used have used the tools of ref. [54]. 

to give the solution Qn 2 = 0.214, shown in Figs. 2-3. We show the results only for the 
scenario II, where C22 = sm2 6 ( see table 1). The predicted cross sections in all cases shown 
are well below the observational limit showing that the spin-independent constraints are 
clearly weaker than the spin-dependent ones despite the coherence enhancement. 

To conclude this section we show the constraints coming from the indirect dark mat- 
ter searches in Super-Kamiokande and IceCube experiments. These limits correspond to 
spin-dependent interactions and are expressed in terms of the standard WIMP-proton inter- 
actions. From Eqn. (4.4) and Table(2) we get the necessary spin-dependent WIMP-proton 
cross section in our model: 

al = dbC p $ = sin 4 ^^[a p (5 p )] 2 J 4 1 . (4.7) 

In our calculations we have used the EMC measurement value for WIMP nucleon spin factor 
ai = 0.46 and of course (S p ) = 0.5 and J = 1/2. In Fig. (5) we show the the predicted 
WIMP-proton cross section in our model together with the constraints as a function of 
the WIMP mass mi . For each value of rri2 the mixing angle sin 9 is again chosen to give 
the correct dark matter density parameter value Oat 2 = 0.214. We have plotted the cross 
sections for mu = 200 GeV and p\2 = +1 and for ms = 500 GeV and p\% = — 1 in the 
scenario II. The other cases shown in Figs. 2-3 are somewhat less constrained. Although 
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m 2 (GeV) m 2 (GeV) 

Figure 5: Shown are spin-dependent WIMP-proton cross sections. Dash dotted blue line shows the 
Super-Kamiokande sensitivity. Red dashed line corresponds to the present IceCube (hard W + W~ 
channel) limit. Black solid lines show the predicted cross sections corresponding to Hn 2 = 0.214 in 
scenario II for m# = 200 GeV and pi2 = +1 (left panel) and for m# = 500 GeV and pi2 = — 1 
(right panel). Again, the standard, no mixing Majorana cross section is shown for comparison with 
the dotted gray line. Finally, the green (light) dashed line shows the predicted IceCube DeepCore 
extension sensitivity [56, 48, 57]. 

the Super-Kamiokande constraint, shown by the blue dash-dotted line, does cut the allowed 
parameter space, the limit is weaker than that coming from the XENON10 experiment. 
The constraint coming from the IceCube experiment, shown by the red dashed (short) line, 
is particularly interesting, because it is much more stringent than the XENON10 limit in 
the high mass region. Yet our model escapes the current IceCube constraint, essentially due 
to the light Higgs final state interactions which force a large suppression on the acceptable 
mixing angles (see section 3). However, with its projected sensitivity after the planned 
DeepCore extension [56, 48], the IceCube experiment will be able to probe a large fraction 
the parameter space available in our model. This is shown by the green dashed (long) line 
in the plot. 

4.3 Oblique constraints 

As we stated already in Sec. 2.4, the required formulas for the evaluation of the oblique 
corrections can be found from [21] and we will not rewrite them here but simply show the 
relevant results. For a fixed phase p±2, we consider only the parameter values (771.2, sin #) 
constrained to satisfy f2jv 2 (7712 , sin 9) ~ 0.214. Then, for any given such pair of values, we 
evaluate the resulting constraints on the remaining mass parameters m\ and rriE when 
the precision parameters S and T are required to lie within limits compatible with present 
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observations. This approach makes sense because, as we have explained above, Q,n 2 depends 
only very weakly on m\ and the- We will use a rather liberal estimate \S\ < 0.2 and 
> T > 0.5. Since the neutrino masses originate from non-renormalizable operators, some 
of the vacuum polarizations required for the calculation of oblique corrections depend on 
the renormalization scale. This scale dependence affects only the polarizations where we 
have neutrinos with masses m a and nib (a,, b = 1,2) in the loops. We have chosen the scale 
as 1.5mi as discussed in [21]. We also require the masses to satisfy ni2 < mi, the- 

In figure 6 we show the constraints from oblique corrections for three different values 
of (777.2, sin 9) which were chosen to probe both small and large mixing and light and heavier 
values of mi. The results corresponding to p\2 = +1 and p\2 = —1 are shown separately in 
left and right panels of the Figure, respectively. From these figures we see that the oblique 
corrections do provide nontrivial constraints on the spectra of the fourth generation leptons 
if the heavy neutrino is required to be responsible for the dark matter abundance of the 
universe. This is so since the oblique corrections are sensitive to all masses in the leptonic 
sector. Moreover there is sensitivity to the mixing angle, because the WIMP N2 contributes 
with a coupling proportional to sin 9 and the heavier mass eigenstate N\ with a coupling 
proportional to ~ cos#. On the other hand we observe that there is still substantial room 
in the parameter space for these masses. In particular the charged lepton can be relatively 
heavy. However, a clear pattern relating the mass of the charged lepton and the heavier 
neutrino eigenstate linearly can be observed. As a rough estimate for our case, niE ~ 2m\ 
allowing for a mild dependence of the slope on the mixing angle 4 . This result is compatible 
with the corresponding results in [21]. Also, the results for positive and negative p\2 are 
qualitatively similar. The hierarchy m,2 < mi,m£ which we have imposed is easy to 
observe from the figures. Also, since the coupling of the charged lepton to electroweak 
currents is proportional to cos 9, we would expect that for small values of sin 9 the oblique 
corrections would favor heavy charged lepton due to large value of the coupling of this 
state; this feature is observed in the above figures and appears more pronounced for the 
case p\2 = —1. Also, due to the interplay of the angles, the most accessible (in view of the 
colliders) mass ranges seem to correspond to sin# ~ 1/2 which "optimizes" the couplings 
of all states to be of similar magnitude. 

5. Conclusions 

We have considered a possibility that the dark matter observed in the universe is a massive 
superweakly interacting Majorana neutrino suggested by recently discovered technicolor 
models for the electroweak symmetry breaking [58, 14]. The spectrum of new particles 
predicted by the model contain a new massive 4th family lepton doublet and a singlet 
sterile state. Our WIMP is the lighter eigenstate in the mixture of the doublet neutrino 
and the sterile state, made stable by a discrete Z2 symmetry under the exchange of any DM- 
sector fields. The mass scale of the new particle spectrum is predicted to be of the order of 
a few hundred GeV by the constraints coming from electroweak precision experiments. The 

4 Note that in section 3 we used rriE = mi. However, we checked that changing the mass hierarchy to 
m,E — 2m 1 has little effect on fijv 2 , as expected. 
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Figure 6: Allowed values of mi and the as constraints \S\ < 0.2 and > T > 0.5 are imposed in ad- 
dition to the constraint £In 2 ( m 2, sin 0). In left panel pi2 = +1 and the three shaded domains corre- 
spond to (sin e,m 2 /m z ) = (0.282,0.55) (lowermost) (0.55,1.29) (top) and (0.223,2.19) (middle). In 
the right panel pi2 = —1 and the three shaded domains correspond to (sin 9, mijvnz) = (0.282, 0.55) 
(lowermost) (0.59,1.29) (top) and (0.207,3.28) (middle). 



relic density of such a neutrino can fall to the observationally determined range £In 2 ~ 0.2 
for a wide range of model parameters, the most important of which are the WIMP and the 
lightest Higgs masses ni2 and run and the active sterile neutrino mixing angle 9. Results 
are also sensitive on the chosen mass generation scenario, of which we have considered two 
distinct possibilities here. In the first scenario all masses are created using only the light 
doublet Higgs composite, while in the second scenario the i?-chiral mass is created by a 
new singlet scalar field. The second scenario allows a much larger range of WIMP masses. 
Yet another parameter in our model is the relative phase P1P2 between the mass eigenstates 
(where pi are the phases necessary to make the masses positive). We pointed out that the 
usual Majorana Feynman rules have to be generalized to compute the necessary matrix 
elements in the presence of such nontrivial mixing phases. The accepted parameter region 
forms a thin surface 9 = 8(rri2, mn , P12), of which we have shown several examples in the 
(m2,sin#)-plane. 

We also computed the constraints on our model parameters coming from the LEP- 
limits on invisible ^-decays, from the corrections to the oblique parameters and from 
direct and indirect dark matter searches. The most stringent constraint on the WIMP 
mass ni2 and mixing angle 9 currently come from the cryogenic XENON10 experiment, 
which already excludes parts of the acceptable parameter space. We also computed the 
indirect bounds coming from Super-Kamiokande and IceCube, finding that the projected 
DeepCore extension of the IceCube experiment will be able to probe a large part of the 
parameter space of our model. However, the most stringent future constraints on the 
model parameters will be set by the upcoming XENON100 and XENON1T experiments. 
Further to these limits, the oblique constraints were shown to be satisfied within narrow 
strips in the plane (mg, mi) of the heavy lepton E and the heavier neutral state N\ masses 
such that the exact position and the area of the strip depends on the given acceptable 
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(m2, sin (9)-pair. Currently it appears that the direct and indirect DM-searches will have 
the better reach to the viability of the model than does the precision electroweak data. In 
particular the forthcoming XENON100 updates will have the sensitivity to rule out almost 
the entire parameter space of our model, or in the more favourable case to detect a WIMP 
with the characteristics proposed by the Minimal Walking Technicolor model. 
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6. Appendix: Cross sections and unitarity 

Computing Majorana cross sections is much more complicated than finding the correspond- 
ing quantities for Dirac particles. Moreover, as shown in section 2.2, the Majorana phases 
enter into the cross sections in quite a nontrivial way with mixing fields. We therefore 
believe that it is worthwhile to present some of the details of our calculations here. While 
most cross sections used in this paper are too long to be written down explicitly, their 
matrix elements can be expressed in a fairly compact form. Even computing these matrix 
elements requires lengthy reductions and combination of many independent contractions, 
and from these results one can already appreciate the nontrivial effects of the mixing phases. 
Matrix elements may also be used as a starting point for other problems involving these 
processes. According to our definition (2.3) N c = pN. This means that the phases pi must 
appear in the field operators, which now become 



where the sum is over the helicities h. This operator gives rise to the following nontrivial 
contractions in the momentum space Feynman rules: 



Only the first two of these will be needed in our computations because in all our processes 
neutrinos appear in the initial state only. In Majorana case one also has four different 
internal contractions, or propagators: 




(6.1) 



a-hpiNi = PiV hi (p) and a hpi Ni 



P*Vhi(p) 
Uhi{p) ■ 



(6.2) 



Ni(x)N i(y) 
Ni(x)Ni(y) 
Ni(x)N t (y) 
Ni(x)Ni(y) 



iSi(x - y) 

iSi{x - y){-piC) 

(pWiSiix-y) 

(pWiSiix-yX-piC), 



(6.3) 
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where iSi{x — y) is the standard Dirac propagator for the mass eigenstate I, and C is 
the usual charge conjugation matrix. The recipe for the cross section calculations is very 
simple. One writes down the S'-matrix element in the standard way, forms all possible 
contractions and uses Eqs. (6.2) and (6.3) to reduce these contractions in terms of Dirac 
matrices. Although a great number of independent contractions are generated in this way, 
the final results of the combinatorics results in rather simple expressions for the matrix 
elements. One then might expect that there should be an effective set of Feynman rules 
that would allow reducing the amount of combinatorics, and several papers have indeed 
put forth such rules, see for example [59]. Unfortunately none of the existing papers offers 
rules that could be used in the case of mixing fields with nontrivial phase factors piPj / 1. 
We believe that simplified rules can be expressed also for mixing fields, but this work will 
be pursued elsewhere [60]. 

For the Gauge boson final states N 2 N 2 — ► CC, where C = W,Z, the matrix elements 
can be generically written in the form 

M CC = ^^1)^2) v N2 ( Pl )P^cT u nM) , (6-4) 

where is the gauge boson polarization vector. The reaction N 2 N 2 — > WW is mediated 
by t- and u-channel -E-lepton exchanges and s-channel Z- and Higgs boson exchanges. In 
each channels one has two independent Majorana contractions. Combining all contractions 
the final matrix element can be compactly written in the form (6.4), with 

= ^ (ASV 1 ^ - Ah" 2 + i#rV 2 (fc - Ah" 1 ) 

+ sin 2 9D z ^ x V^_ ku _ k2 - 2m 2 C\ 2 D H g^ , (6.5) 
where the three-gauge-boson coupling factor is 

V k ]^_ ku _ k2 = (k 2 - h) x g^ 2 ~ (h + 2fe)^ + (k 2 + 2hrg^ (6.6) 
and the propagator factors Dx are given by (3.7) for X = Z,H and 



1 

a + 



for a = t,u. Like fermion final state processes, the WW-process do not involve neutrino 
mixing vertices, and so the result involves no nontrivial phases and it agrees with the one 
obtained by use of the standard Majorana Feynman rules. Obviously, the expression for 
the total cross section computed from Mww must be a very lengthy, and we shall not 
present it here. (For an analogous expression of an annihilation of Dirac neutrinos see 
ref. [30].) Let us however consider explicitly the limit s^oo. The leading term of the total 
cross section now becomes: 



G 2 m 2 (\ 
<tn 2 n 2 ^w+w- (s > m x ) ~ -^(C^ - sin 2 6) 2 + O ( - log . 



(6.8) 
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Obviously this cross section is unitary in the scenario II but not in scanario I. The unitarity 
of scenario II is actually an exception and nonunitarity is a common feature for all our gauge 
and Higgs boson final states with mixing Majorana fields. However, the nonunitary is not 
a problem here, because our low-energy theory giving rise to the masses is only an effective 
one, known to be broken at the ~TeV region. 

In the process N2N2 — > ZZ one has N\- and ^-mediated mediated t- and u-channel 
processes an iJ-mediated s-channel processes. Here each t— and u— channel process con- 
sists of four independent contractions. Including two nontrivial s-channel contractions, one 
has to work out and combine 18 independent matrix elements. The final result can again 
be written in the form (6.4), with 

= 2cosW < Sin4 9 (A ~ h + ™ 2)yi2 + ^ (A ~ h + m2h ^ 

+ sin 2 flcos 2 6> [AVHA ~h +m 1 h>* 2 + J D 2 7 ' l2 (^ 1 - jt 2 + m 1 )^ 1 ] 

+ 8m 2 C^D H g^y (6.9) 

where fh\ = pv2"m\. Depending on the sign of the relative phase p\2 = p\P2 the effective 
mass of the heavier field N 2 can be either positive or negative in this expression. Different 
choices for this phase lead to physically different cross sections and eventually very different 
predictions for the relic density. To our knowledge this is the first time when Majorana 
phases have been observed to have such a large effect on physical observables. In the limit 
s — > 00 the cross section for process N2N 2 ^ZZ in the scenario II becomes: 

Q2 

<jn 2 n 2 ->zz(s > rax) « — - sin 4 9 cos 4 9 (m| - pi2m\) 2 

= 91 ^ (610 ) 

8vr (M L - M R ) 2 + 4m 2 D ' K ' ' 

Again the nonunitary behaviour is manifest. However, unitarity is again restored in the 
limit of no mixing, or equivalently, for a vanishing Dirac mass term. This restoration of 
unitarity at zero mixing is a generic result for all cross sections in the scenario II, following 
from the fact that in this limit the WIMP becomes a purely sterile state, which is con- 
nected to the dangerous gauge boson degrees of freedom only through the Dirac mass term. 

The annhilation to higgs bosons, N 2 N2 — > HH processes through both N\ and N2- 
mediated t- and u-channel diagrams as well as an s-channel Higgs exchange and the 
ft./iiV2-/V2-contact interaction. After some algebra the matrix element can be written as 



Mhh = -^-^rV N2 {pi)piA HH u N2 {p2) , (6.11) 



where 



A HH = 2{C%2? [D 2 t {2m 2 - jt 2 ) + D 2 u (2m 2 - JtJ] 

+ U°i2) 2 [Dl(m 2 + mi - h) + D l( m 2 + mi - h)] 

+ 3C^D H + —C£. (6.12) 
m 2 rri2 
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where again rh\ = pvimx. The coefficients C% 2 , C\ 2 an d C\2 are given in the Table 1 and 
the propagators D a in Eq. (6.7). 

Finally, the annihilations to the mixed final state ZH are mediated by the t- and 
it-channel N\ and N 2 exchanges, as well as by an s-channel Z-exchange diagram. We find: 

• o 

MzH = 2 cos 2 Bw Wz^^ ^ N ^ Pl) P iA hz u n 2 {P2) , (6.13) 

where 

A% H = sin 2 0C* [Dh^lHh ~h + m2) + D 2 u {i> 2 -jt 1 + m 2 ) 7 V] 

+ \ sin2^ 12 Cf 2 [AVT 5 ^ - h + mi) + £>i(& - h + ™ihV] 
2M 2 

+ ^ sin 2 0D Z ^-/ 5 ■ (6.14) 

m 2 

When matrix elements are given in these compact form, it is a straightforward to compute 
the cross sections using some algebraic manipulation program. We reduced all the traces 
using FEYNCALC, after which the resulting cross sections were transported to a dedicated 
fortran code which evaluated (va) at a given temperature T from the integral expression 
(3.4). 
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